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Abstract. We review our recent development of family replicated gauge group model, 
which generates the Large Mixing Angle MSW solution. The model is based on each 
' family of quarks and leptons having its own set of gauge fields, each containing a 

replica of the Standard Model gauge fields plus a (B — L)-coupled gauge field. A fit 
of all the seventeen quark-lepton mass and mixing angle observables, using just six 
new Higgs field vacuum expectation values, agrees with the experimental data order of 
magnitudewise. However, this model can not predict the baryogenesis in right order, 
therefore, we discuss further modification of our model and present a preliminary result 
of baryon number to entropy ratio. 
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1 Introduction 



We have previously attempted to fit all the fermion masses and their mixing 
pi 1 angles including baryogenesis (|] in a model without supersymmetry or 

grand unification. This model has the maximum number of gauge fields con- 
sistent with maintaining the irreduciblity of the usual Standard Model fermion 
<~] \ representations, added three right-handed neutrinos. The predictions of this pre- 

vious model are in order of magnitude agreement with all existing experimental 
data, however, only provided we use the Small Mixing Angle MSW [Q (SMA- 
MSW) solution. But, for the reasons given below, the SMA-MSW solution is 
now disfavoured by experiments. So here we review a modified version of the 
previous model, which manages to accommodate the Large Mixing Angle MSW 
(LMA-MSW) solution for solar neutrino oscillations using 6 additional Higgs 
fields (relative to the Standard Model) vacuum expectation values (VEVs) as 
adjustable parameters. 

A neutrino oscillation solution to the solar neutrino problem and a favouring 
of the LMA-MSW solution is supported by SNO results ||: The measurement 
of the 8 B and hep solar neutrino fluxes shows no significant energy dependence 
of the electron neutrino survival probability in the Super-Kamiokande and SNO 
energy ranges. 

Moreover, the important result which also supports LMA-MSW solution on 
the solar neutrino problem, reported by the Super-Kamiokande collaboration Q , 
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that the day-night asymmetry data disfavour the SMA-MSW solution at the 95% 
C.L.. 

In fact, global analyses @,|]||[l(| of all solar neutrino data have confirmed 
that the LMA-MSW solution gives the best fit to the data and that the SMA- 
MSW solution is very strongly disfavoured and only acceptable at the 3<r level. 
Typical best fit values of the mass squared difference and mixing angle param- 
eters in the two flavour LMA-MSW solution are Am.Q w 4.5 x 10~ 5 eV 2 and 
tan 2 9 Q w 0.35. 

This paper is organised as follows: In the next section, we present our gauge 
group - the family replicated gauge group - and the quantum numbers of fermion 
and Higgs fields. Then, in section 3 we discuss our philosophy of all gauge- and 
Yukawa couplings at Planck scale being of order unity. In section 4 we address 
how the family replicated gauge group breaks down to Standard Model gauge 
group, and we add a small review of see-saw mechanism. The mass matrices 
of all sectors are presented in section 5, the renormalisation group equations - 
renormalisable and also 5 dimensional non-renormalisable ones - are shown in 
section 6. The calculation is described in section 7 and the results are presented in 
section 8. We discuss further modification of our model and present a preliminary 
results of baryon number to entropy ratio in section 9. Finally, section 10 contains 
our conclusion. 

2 Quantum numbers of model 

Our model has, as its back-bone, the property that there are generations (or 
families) not only for fcrmions but also for the gauge bosons, i.e., we have a 
generation (family) replicated gauge group namely 

^ASMG., x U(1) B _ LA ) , (1) 

where SMG denotes the Standard Model gauge group = SU(3) x SU(2) x [7(1), 
x denotes the Cartesian product and i runs through the generations, i — 1,2,3. 

Note that this family replicated gauge group, eq. ([!]), is the maximal gauge 
group under the following assumptions: 

• It should only contain transformations which change the known 45 (= 3 
generations of 15 Weyl particles each) Weyl fermions of the Standard Model 
and the additional three heavy see-saw (right-handed) neutrinos. That is our 
gauge group is assumed to be a subgroup of t/(48). 

• We avoid any new gauge transformation that would transform a Weyl state 
from one irreducible representation of the Standard Model group into an- 
other irreducible representation: there is no gauge coupling unification. 

• The gauge group does not contain any anomalies in the gauge symmetry - 
neither gauge nor mixed anomalies even without using the Green-Schwarz 
anomaly cancelation mechanism. 

• It should be as big as possible under the foregoing assumptions. 
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Table 1. All U(l) quantum charges in the family replicated model. The symbols for 
the fermions shall be considered to mean "proto" -particles. Non-abelian representations 
are given by a rule from the abelian ones (see Eq. (H)). 
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The quantum numbers of the particles/fields in our model are found in table 
1 and use of the following procedure: In table 1 one finds the charges under the 
six U(l) groups in the gauge group |[ Then for each particle one should take 
the representation under the SU(2)i and SU(3)i groups (i = 1, 2, 3) with lowest 
dimension matching to yt/2 according to the requirement 

| + | + | = (modi) , (2) 

where U and di are the triality and duality for the i'th proto-generation gauge 
groups SU(3)i and SU(2)i respectively. 

3 The philosophy of all couplings being order unity 

Any realistic model and at least certainly our model tends to get far more fun- 
damental couplings than we have parameters in the Standard Model and thus 
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pieces of data to fit. This is especially so for our model based on many U(l) 
charges because we take it to have practically any not mass protected par- 
ticles one may propose at the fundamental mass scale, taken to be the Planck 
mass. Especially we assume the existence of Dirac fermions with order of fun- 
damental scale masses needed to allow the quark and lepton Weyl particles to 
take up successively gauge charges from the Higgs fields VEVs. So unless we 
make assumptions about the many coupling constants and fundamental masses 
we have no chance to predict anything. Almost the only chance of making an 
assumption about all these couplings, which is not very model dependent, is to 
assume that they are all of order unity in the fundamental unit. This is the 
same type of assumption that is really behind use of dimensional arguments to 
estimate sizes of quantities. A procedure very often used successfully. If we re- 
ally assumed every coupling and mass of order unity we would get the effective 
Yukawa couplings of the quarks and leptons to the Weinberg-Salam Higgs field 
to be also of order unity what is phenomenologically not true. To avoid this 
prediction we then blame the smallncss of all but the top- Yukawa coupling on 
smallness in fundamental Higgs VEVs. That is to say we assume that the VEVs 
of the Higgs fields in Table 1, p, ui, T, W, x, 4»b-l and 4*ws are (possibly) very 
small compared to the fundamental/Planck unit, and these are the quantities 
we have to fit. 

Technically we implement these unknown - but of order unity according to 
our assumption - couplings and masses by a Monte Carlo technique: we put 
them equal to random numbers with a distribution dominated by numbers of 
order unity and then perform the calculation of the observable quantities such 
as quark or lepton masses and mixing angles again and again. At the end we 
average the logarithmic of these quantities and exponentiate them. In this way 
we expect to get the typical order of magnitude predicted for the observable 
quantities. In praxis we do not have to put random numbers in for all the many 
couplings in the fundamental model, but can instead just provide each mass 
matrix element with a single random number factors. 

After all a product of several of order unity factors is just an order unity 
factor again. To resume our model philosophy: Only Higgs field VEVs are not 
of order unity. We must be satisfied with order of magnitude results. 

4 Breaking of the family replicated gauge group to the 
Standard Model 

The family replicated gauge group broken down to its diagonal subgroup at 
scales about one or half order of magnitude under the Planck scale by Higgs 
fields - W, T, to, p and x (in Table 0): 

i= * 2i3 (SMG l x U{1) B _ LA ) -> SMG x U(1) B -l ■ (3) 

This diagonal subgroup is further broken down by yet two more Higgs fields 
- the Weinberg-Salam Higgs field 4>ws and another Higgs field 4>b-l — to 
SU{2) x U(l) em . 
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4.1 See-saw mechanism 

See-saw mechanics is build into our model to fit the scale of the neutrino os- 
cillations, i.e., we use the right-handed neutrinos with heavy Majorana masses 
(10 11 GeV). 

In order to mass-protect the right-handed neutrino from getting Planck scale 
masses, we have to introduce <p B -L which breaks the B — L quantum charge 
spontaneously, and using this new Higgs filed we are able to deal the neutrino 
oscillations, i.e., to fit the scale of the see-saw particle masses. However, due 
to mass-protection by the Standard Model gauge symmetry, the left-handed 
Majorana mass terms should be negligible in our model. Then, naturally, the 
light neutrino mass matrix - effective left-left transition Majorana mass matrix 
- can be obtained via the see-saw mechanism j[2| : 

M eB ttM? M^iMff . (4) 



5 Mass matrices 



Using the U(l) fermion quantum charges and Higgs field (presented in Table [[]) 
we can calculate the degrees of suppressions of the left-right transition - Dirac 
mass - matrices and also Majorana mass matrix (right-right transition). 

Note that the random complex order of unity numbers which are supposed to 
multiply all the mass matrix elements are not represented in following matrices: 
the up-type quarks: 
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and the Majorana (right-handed) neutrinos: 
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(5) 



(6) 
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6 Renormalisation group equations from Planck scale to 
week scale via see-saw scale 

It should be kept in mind that the effective Yukawa couplings for the Weinberg- 
Salam Higgs field, which are given by the Higgs field factors in the above mass 
matrices multiplied by order unity factors, are the running Yukawa couplings 
at a scale near the Planck scale. In this way, we had to use the renormalisation 
group (one-loop) /3-functions to run these couplings down to the experimentally 
observable scale which we took for the charged fermion masses to be compared to 
"measurements" at the scale of 1 GeV, except for the top quark mass prediction. 
We define the top quark pole mass: 

M t = mt {M)(l + ^T) , (10) 

where we put M = 180 GeV for simplicity. 

We use the one-loop (3 functions for the gauge couplings and the charged 
fermion Yukawa matrices [|l3) as follows: 

dt 10 J1 

, R 2 d 92 19 3 

1d7t — — = q, 

dt 16 y2 

1C 2^53 „ 3 

™?^£ = \ -Y D <Y B )t)Y v 



16n 2 ^ = 3 ~(Y D (Y n ?-Y P (Y v )i)Y D 

Ys ~ (i-9? + ^2+8.g 3 2 )} Y D 

y s =Tr(3y [ ty c , + 3Y D ty D + y E ty E ) , 

where t = ln/x. 

By calculation we use the following initial values of gauge coupling constants: 

U(l): 9l (M z ) = 0.462 , 5 i(M Planck ) = 0.614 (12) 
5(7(2): g 2 (M z ) = 0.651 , g 2 (M Planck ) = 0.504 (13) 
5(7(3) : g 3 (M z ) = 1.22 , 53 (M Planck ) = 0.491 (14) 
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6.1 The renormalisation group equations for the effective neutrino 
mass matrix 

The effective light neutrino masses are given by an irrelevant, nonrenormal- 
isable (5 dimensional term) - effective mass matrix M e ff - for which the running 
formula is Gj : 

16n2 ^f L = (- 3 52+2A + 2r s ) M eff - | (M eff (Y E Y^ f + (Y B Yj) M off ) , (15) 

where A is the Weinberg-Salam Higgs self-coupling constant and the mass of the 
Standard Model Higgs boson is given by Mfj — A {4>ws) 2 ■ We just for simplicity 
take Mh = 115 GeV thereby we ignore the running of the Higgs self-coupling 
and fixed as A = 0.2185. 

Note that the renormalisation group equations are used to evolve the effective 
neutrino mass matrix from the see-saw sale, set by {4>b~l) in our model, to 
1 GeV. 

7 Method of numerical computation 

In the philosophy of order unity numbers spelled out in section || we evaluate the 
product of mass-protecting Higgs VEVs required for each mass matrix element 
and provide it with a random complex number, Ajj, of order one as a factor 
taken to have Gaussian distribution with mean value zero. But we hope the 
exact form of distribution does not matter much provided we have (In = 0. 
In this way, we simulate a long chain of fundamental Yukawa couplings and 
propagators making the transition corresponding to an effective Yukawa coupling 
in the Standard Model and the parameters in neutrino sector. In the numerical 
computation we then calculate the masses and mixing angles time after time, 
using different sets of random numbers and, in the end, we take the logarithmic 
average of the calculated quantities according to the following formula: 

(m) = exp ^ . (16) 

Here (to) is what we take to be the prediction for one of the masses or mixing 
angles, m, is the result of the calculation done with one set of random number 
combinations and N is the total number of random number combinations used. 

Since we only expect to make order of magnitude fits, we should of course 
not use ordinary \ 2 defined form the experimental uncertainties by rather the 
X 2 that would correspond to a relative uncertainly - an uncertain factor of order 
unity. Since the normalisation of such a x 2 is n °t so easy to choose exactly we 
define instead a quantity which we call the goodness of fit (g.o.f). Since our 
model can only make predictions order of magnitudewise, this quantity g.o.f. 
should only depend on the ratios of the fitted masses and mixing angles to the 
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experimentally determined masses and mixing angles: 

where m cxp are the corresponding experimental values presented in Table ||. 

We should emphasise that we do not adjust the order of one numbers by 
selection, i.e., the complex random numbers are needed for only calculational 
purposes. That means that we have only six adjustable parameters - VEVs 
of Higgs fields - and, on the other hand, that the averages of the predicted 
quantities, (m), are just results of integration over the "dummy" variables - 
random numbers - therefore, the random numbers are not at all parameters! 

Strictly speaking, however, one could consider the choice of the distribution 
of the random order unity numbers as parameters. But we hope that provided 
we impose on the distribution the conditions that the average be zero and the 
average of the logarithm of the numerical value be zero, too, any reasonably 
smooth distribution would give similar results for the (m) values at the end. In 
our early work we did see that a couple of different proposals did not make 
too much difference. 



(17) 



8 Results 

We averaged over N = 10, 000 complex order unity random number combina- 
tions. These complex numbers are chosen to be a number picked from a Gaussian 
distribution, with mean value zero and standard deviation one, multiplied by a 
random phase factor. We put them as factors into the mass matrices (eqs. 
Then we computed averages according to eq. (|l^) and used eq. ( |l7| ) as a \ 2 t° 
fit the 6 free parameters and found: 

(4> ws ) = 246 GeV, (4> B _ L ) = 1.64 x 10 11 GeV, (u) = 0.233, 

(p) = 0.246 , (W) = 0.134 , (T) = 0.0758 , ( X ) = 0.0737 , (18) 

where, except for the Weinberg-Salam Higgs field and (4>b-l), the VEVs are 
expressed in Planck units. Hereby we have considered that the Weinberg-Salam 
Higgs field VEV is already fixed by the Fermi constant. The results of the best 
fit, with the VEVs in eq. (]l8|), are shown in Table ^ and the fit has g.o.f. = 3.63. 

We have 11 = 17 — 6 degrees of freedom - predictions - leaving each of 
them with a logarithmic error of v/3-63/11 ~ 0.57, which is very close to the 
theoretically expected value 64% nM. This means that we can fit all quantities 
within a factor exp 63/11^ ~ 1.78 of the experimental values. 

From the table |2j the experimental mass values are a factor two higher than 
predicted for down, charm and for the Cabibbo angle V us while they are smaller 
by a factor for strange and electron. Thinking only on the angles and masses 
(not squared) the agreement is in other cases better than a factor two. 

Experimental results for the values of neutrino mixing angles are often pre- 
sented in terms of the function sin 2 29 rather than tan 2 9 (which, contrary to 



Family replicated fit of all quark and lepton masses and mixings 



9 



Table 2. Best fit to conventional experimental data. All masses are running masses 
at 1 GeV except the top quark mass which is the pole mass. Note that we use the 
square roots of the neutrino data in this Table, as the fitted neutrino mass and mixing 
parameters (to), in our goodness of fit (g.o.f.) definition, eq. (P-7|). 
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sin 2 20, does not have a maximum at 6 — n/4 and thus still varies in this re- 
gion). Transforming from tan 2 9 variables to sin 2 28 variables, our predictions for 
the neutrino mixing angles become: 

sin 2 26» =0.61, (19) 
sin 2 20 atm = 0.96, (20) 
sin 2 26»i 3 =0.17. (21) 

We also give here our predicted hierarchical neutrino mass spectrum: 

(22) 
(23) 
(24) 

Our agreement with experiment is excellent: all of our order of magnitude 
neutrino predictions lie inside the 99% C.L. border determined from phenomeno- 
logical fits to the neutrino data, even including the CHOOZ upper bound. Our 
prediction of the solar mass squared difference is about a factor of 2 larger than 
the global data fit even though the prediction is inside of the LMA-MSW re- 
gion, giving a contribution to our goodness of fit of g.o.f. « 0.14. Our CHOOZ 
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angle also turns out to be about a factor of 2 larger than the experimental limit 
at 90% C.L., delivering another contribution of g.o.f. ~ 0.14. In summary our 
predictions for the neutrino sector agree extremely well with the data, giving a 
contribution of only 0.34 to g.o.f. while the charged fermion sector contributes 
3.29 to g.o.f.. 

8.1 CP violation 

Since we have taken our random couplings to be - whenever allowed - complex 
we have order of unity or essentially maximal CP-violation so a unitary triangle 
with angles of order one is a success of our model. After our fitting of masses 
and of mixings we can simply predict order of magnitudewise of CP-violation 
in e.g. K° — K° decay or in CKM and MNS mixing matrices in general. 

The Jarlskog area J CP provides a measure of the amount of CP violation 
in the quark sector ]Tr| and, in the approximation of setting cosines of mixing 
angles to unity, is just twice the area of the unitarity triangle: 

Jcp = V us V c b V ub sin 5 , (25) 

where S is the CP violation phase in the CKM matrix. In our model the quark 
mass matrix elements have random phases, so we expect S (and also the three 
angles a, (5 and 7 of the unitarity triangle) to be of order unity and, taking an 
average value of | sin <5| f=a 1/2, the area of the unitarity triangle becomes 

Jcp « \ Vus Vcb V ub ■ (26) 

Using the best fit values for the CKM elements from Table ||, we predict J CP « 

3.1 x 10~ 6 to be compared with the experimental value (2 — 3.5) x 10~ 5 . Since 
our result for the Jarlskog area is the product of four quantities, we do not 
expect the usual ±64% logarithmic uncertainty but rather ±a/4 • 64% = 128% 
logarithmic uncertainty. This means our result deviates from the experimental 
value by ln( g [^°Ig )/1.28 = 1.7 "standard deviations". 

The Jarlskog area has been calculated from the best fit parameters in Table 
2, it is also possible to calculate them directly while making the fit. So we have 
calculated J CP for N = 10, 000 complex order unity random number combina- 
tions. Then we took the logarithmic average of these 10, 000 samples of J CP and 
obtained the following result: 

J OP = 3.1 x 10~ 6 , (27) 
in good agreement with the values given above. 

8.2 Neutrinoless double beta decay 

Another prediction, which can also be made from this model, is the electron 
"effective Majorana mass" - the parameter in neutrinoless beta decay - defined 
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by: 
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(28) 



where rm are the masses of the neutrinos v% and U e i are the MNS mixing matrix 
elements for the electron flavour to the mass eigenstates i. We can substitute 
values for the neutrino masses from eqs. ( p2] - |2^ ) and for the fitted neu- 
trino mixing angles from Table || into the left hand side of eq. (^8|) . As already 
mentioned, the CP violating phases in the MNS mixing matrix are essentially 
random in our model. So we combine the three terms in eq. (|2^) by taking the 
square root of the sum of the modulus squared of each term, which gives our 
prediction: 

| (m) | « 3.1 x 1CT 3 eV. (29) 

In the same way as being calculated the Jarlskog area we can compute using 
N = 10, 000 complex order unity random number combinations to get the |(to)|. 
Then we took the logarithmic average of these 10, 000 samples of | (to) | as usual: 



\{m}\ = 4.4 x 10~ 3 eV. 

This result does not agree with the central value of recent result 
from the Heidelberg-Moscow collaboration Wn. 



(30) 



"evidence" 



9 Baryogenesis via Lepton Number Violation 

Having now a well fitted model giving orders of magnitude for all the Yukawa 
couplings and having the see-saw mechanism, it is obvious that wc ought to 
calculate the amount of baryons Yb relative to entropy being produced via the 
Fukugita-Yanagida mechanism Jl9]. According to this mechanism the decay of 
the right-handed neutrinos by CP-violating couplings lead to an excess of the 
B — L charge (meaning baryon number minus lepton number), the relative ex- 
cess in the decay from Majorana neutrino generation number i being called 6j. 
This excess is then immediately - and continuously back and forth - being con- 
verted partially to a baryon number excess, although it starts out as being a 
lepton number L asymmetry, since the right-handed neutrinos decay to leptons 
and Weinberg-Salam Higgs particles. It is a complicated discussion to estimate 
to what extend the B — L asymmetry is washed out later in the cosmological 
development, but our estimates goes that there is not enough baryon number 
excess left to fit the Big Bang development at the stage of formation of the light 
elements primordially (nuclearsynthesis). 

Recently we have, however, developed a modified version |]l9f of our model - 
only deviating in the right-handed sector - characterized by changing the quan- 
tum numbers assumed for the see-saw scale producing Higgs field 4>b-l m such 
a way that the biggest matrix elements in the right-handed mass matrix (eq. ^ 
becomes the pair of - because of the symmetry - identical off diagonal elements 
(row, column)=(2,3) and (3,2). Thereby we obtain two almost degenerate right- 
handed neutrinos and that helps for making the B — L asymmetry in the decay 
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bigger. In this modified model that turns out to fit the rest of our predictions 
approximately equally well or even better we then get a very satisfactory baryon 
number relative to entropy prediction 

Y B w 2.5 x 1CT 11 . (31) 

In the same time as making this modification of the 4>b-l quantum numbers 
we also made some improvements in the calculation by taking into account the 
running of the Dirac neutrino Yukawa couplings from the Planck scale to the 
corresponding right-handed neutrino scales. Also, we calculated more accurate 
dilution factors than previous our work || . However, foregoing work was based 
on the mass matrices which predicted the SMA-MSW, so we must investigate 
the baryogenesis using the present mass matrices, of course, with the modified 
right-handed Majorana mass matrix. 

10 Conclusion 

We have reviewed our model which is able to predict the experimentally favored 
LMA-MSW solution rather than the SMA-MSW solution for solar neutrino os- 
cillations after careful choice of the U(l) charges for the Higgs fields causing 
transitions between 1st and 2nd generations. However, the fits of charged lepton 
quantities become worse compare to our "old" model that can predict SMA- 
MSW solar neutrino solution. On the other hand, we now can fit the neutrino 
quantities very well: the price paid for the greatly improved neutrino mass matrix 
fit - the neutrino parameters now contribute only very little to the g.o.f. - is a 
slight deterioration in the fit to the charged fermion mass matrices. In particular 
the predicted values of the quark masses to^ and m c and the Cabibbo angle V us 
are reduced compared to our previous fits. However the overall fit agrees with the 
seventeen measured quark-lepton mass and mixing angle parameters in Table 2 
within the theoretically expected uncertainty |lfj] of about 64%; it is a perfect 
fit order of magnitudewisc. It should be remarked that our model provides an 
order of magnitude fit /understanding of all the effective Yukawa couplings of 
the Standard Model and the neutrino oscillation parameters in terms of only 6 
parameters - the Higgs field vacuum expectation values! 
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